Gapless Line for the Anisotropic Heisenberg Spin- 1/2 Chain in a Magnetic Field 

and the Quantum ANNNI Chain 



m 
O 
O 
(N 

Of 
< 
OO 



CO 

c3 



O 

o 



> 

(N 

O 
m 
O 

03 



c 

o 
o 



X 
S3 



Amit Dutta 1 * and Diptiman Sen 2 
1 Institut fur Theoretische Physik und Astrophysik, Universitdt Wurzburg, Am Hubland, 97074 Wiirzburg, Germany 
2 Centre for Theoretical Studies, Indian Institute of Science, Bangalore 560012, India 

(February 2, 2008) 

We study the anisotropic Heisenberg (XYZ) spin- 1/2 chain placed in a magnetic field pointing 
along the x-axis. We use bosonization and a renormalization group analysis to show that the model 
has a non-trivial fixed point at a certain value of the XY anisotropy a and the magnetic field h. 
Hence, there is a line of critical points in the (a, h) plane on which the system is gapless, even though 
the Hamiltonian has no continuous symmetry. The quantum critical line corresponds to a spin-flop 
transition; it separates two gapped phases in one of which the Z2 symmetry of the Hamiltonian is 
broken. Our study has a bearing on one of the transitions of the axial next-nearest neighbor Ising 
(ANNNI) chain in a transverse magnetic field. We also discuss the properties of the model when 
the magnetic field is increased further, in particular, the disorder line on which the ground state is 
a direct product of single spin states. 



PACS number: 75.10.Jm, 64.70.Rh 



I. INTRODUCTION 

One-dimensional quantum spin systems have been 
studied extensively ever since the problem of the 
isotropic Heisenberg spin-1/2 chain was solved exactly by 
Bethe. Baxter later used the Bethe ansatz to solve the 
anisotropic Heisenberg (XYZ) spin-1/2 chain in the ab- 
sence of a magnetic field [1]; the problem has not been an- 
alytically solved in the presence of a magnetic field. Ex- 
perimentally, quantum spin chains and ladders are known 
to exhibit a wide range of unusual properties, including 
both gapless phases with a power-law decay of the two- 
spin correlations and gapped phases with an exponential 
decay [2,3]. There are also two-dimensional classical sta- 
tistical mechanics systems (such as the axial next-nearest 
neighbor Ising (ANNNI) model) whose finite temperature 
properties can be understood by studying an equivalent 
quantum spin-1/2 chain in a magnetic field. The ANNNI 
model has been studied by several techniques, and it was 
believed for a long time to have a floating phase of finite 
width in which the system is gapless [4] . 

Amongst the powerful analytical methods now avail- 
able for studying quantum spin-1/2 chains is the tech- 
nique of bosonization [2,5]. Recently, the XX Z chain in 
a transverse magnetic field [6] and the quantum ANNNI 
model [7] have been studied using bosonization. In this 
paper, we will study the anisotropic XYZ model in a 
magnetic field pointing along the ir-axis. For small values 
of the XY anisotropy a and the magnetic field h, we will 
show in Sec. II that there is a non-trivial fixed point (FP) 
of the renormalization group (RG) in the (a, h) plane; 
the system is gapless on a quantum critical line of points 
which flow to this FP. In Sec. Ill, one of the transitions 
of the ANNNI model will be shown to be a special case of 



our results in which the zz coupling is equal to zero. Our 
results are complimentary to the earlier studies of the 
ANNNI model which indicate a gapless phase of finite 
width. We will present the complete zero temperature 
phase diagram of the ANNNI model which has both a 
gapless phase of finite width as well as a gapless line. 

The gapless line is somewhat unusual because the XY 
anisotropy and the magnetic field both break the contin- 
uous symmetry of rotations in the x — y plane. In Sec. 
IV, we will provide a physical understanding of the gap- 
less line by going to the classical (large S) limit of the 
model; this helps us to identify it as a spin-flop transition 
line. In Sec. V, we will discuss a disorder line which lies 
at a larger value of the magnetic field. In Sec. VI, we 
will briefly comment on the Ising transition which occurs 
at an even larger value of the magnetic field. 



II. BOSONIZATION AND RENORMALIZATION 
GROUP ANALYSIS 

We consider the anisotropic Hamiltonian defined on a 
chain of sites 



H = 



E 



(1 + a) SISl, 



(1 - a) sy n s v 



n+1 



h SI 



(1) 



where the S" are spin-1/2 operators. We will assume 
that the XY anisotropy a and the zz coupling A satisfy 
— 1 < a, A < 1. We can assume without loss of gen- 
erality that the magnitude of the zz coupling is smaller 
than the yy coupling (i.e., |A| < 1 — a), and that the 
magnetic field strength h > 0. The Hamiltonian in 
Eq. (1) is invariant under the global Z2 transformation 
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For a = h = 0, the model is symmetric under rotations 
in the x — y plane and is gapless. The low-energy and 
long-wavelength modes of the system are then described 
by the bosonic Hamiltonian [2,5] 



dx [ (d x ef + {d x 4>f 



(2) 



where v is the velocity of the low-energy excitations 
(which have the dispersion u> = v\k\); v is a function of A. 
(The continuous space variable x and the site label n are 
related through x = nd, where d is the lattice spacing.) 
The bosonic theory contains another parameter called K 
which is related to A by [2,5] 



K 



ir + 2 sin - (A) 



(3) 



K takes the values 1 and 1/2 for A = (which de- 
scribes noninteracting spinless fermions) and A = 1 
(the isotropic antiferromagnet) respectively; as A — > — 1, 
K — ► co. We thus have 1/2 < K < oo. 

In terms of the fields <j> and 9 introduced in Eq. (2), 
the spin operators can be written as [6] 



S n = \Ik d ^ + (-l)"cicos(2V^>) , 



S x n = [c 2 cos(2V^) + (-l) n cs] cos^-6), (4) 

where the Cj are constants given in Ref. 8. The XY 
anisotropy term is given by 



^n^n+l ~ c 4 COs(2 



where C4 is another constant. 

For convenience, let us define the three operators 



(5) 



O x = cos(2Vt^</>)cos( 



e> 2 = cos(2^/— 6») , and 3 = coa(WirK4>) . (6) 

Their scaling dimensions are given by K + 1/AK, 1/K 
and AK respectively. Using Eqs. (4-5), the terms corre- 
sponding to a and h in Eq. (1) can be written as 



H a + H h 



= I dx[ 



(7) 



where we have dropped rapidly varying terms propor- 
tional to (—1)™ since they will average to zero in the 
continuum limit. (We will henceforth absorb the factors 
c 4 (c 2 ) in the definitions of a (h).) We will now study 
how the parameters a and h flow under RG. 

The operators in Eqs. (6) are related to each other 
through the operator product expansion; the RG equa- 
tions for their coefficients will therefore be coupled to 



each other [9]. In our model, this can be derived as fol- 
lows. Given two operators A\ = exp(icti<f) + i(3\6) and 
A 2 = exp(ia20 + ifcO), we write the fields 4> and 9 as 
the sum of slow fields (with wave numbers \k\ < Ae~ dl ) 
and fast fields (with wave numbers Ae~ dl < \k\ < A), 
where A is the momentum cut-off of the theory, and dl is 
the change in the logarithm of the length scale. Integrat- 
ing out the fast fields shows that the product of A\ and 
Ai at the same space-time point gives a third operator 



A, 



o i(ai+a 2 )0+i(/3i+/32) 



with a prcfactor which can be 



schematically written as 



AxA 2 



e -{a 1 a 2 +f3 1 f3 2 )dl/2ir a 



(8) 



If denote the coefficients of the operators Ai in an 
effective Hamiltonian, then the RG expression for d\% / dl 
will contain the term (aia 2 + /3i/32)AiA 2 /27r. Using this, 
we find that if the three operators in Eqs. (6) have co- 
efficients h, a and b respectively, then the RG equations 
are 



dh 
~dl 

da 

df 
db 

~dl 

dJK 
~df 



(2 - K 



AK' 



1 

K 



ah - AKbh 



(2-4*06 + {2K-^ R )h 2 , 



(9) 



where we have absorbed some factors involving v in the 
variables a, b and h. (We will ignore the RG equation for 
v here.) It will turn out that K renormalizes very little in 
the regime of RG flows that we will be concerned with. 
Eqs. (9) have appeared earlier in the context of some 
other problems [10,11]. However, the last two terms in 
the expression for dh/dl were not presented in Ref. 10; 
these two terms turn out to be crucial for what follows. 
Note that Eqs. (9) are invariant under the duality trans- 
formation K <-» 1/AK and a <-» b. 

Let us now consider the fixed points of Eqs. (9). For 
any value of K = K* , a trivial FP is (a*,b*,h*) = 
(0,0,0). Remarkably, it turns out that there is a non- 
trivial FP for any value of K* lying in the range 1/2 < 
K* < l + y/3/2; we will henceforth restrict our attention 
to this range of values. (The upper bound on K* comes 
from the condition 2 - K* - 1 j AK* > 0.) The non- trivial 
FP is given by 



h* = 



V / 2K*(2-K* - 1/AK*) 
2K* + 1 



a* = (K* + l)h* 2 , and b* = 

\ 2 > 2K* 



(10) 



The system is gapless at this FP as well as at all points 
which flow to this FP. One might object that Eqs. (9) 
can only be trusted if a, b and h are not too large, oth- 
erwise one should go to higher orders. We note that the 
FP approaches the origin as K* — > 1 + V3/2 ~ 1.866; 
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from Eq. (3), this corresponds to the zz coupling A = 
— sin[7r(v / 3— 3/2)] ~ —0.666. Thus the RG equations can 
certainly be trusted for K* close to 1.866. For K* = 1, 
the FP is at (a*,b*,h*) = (1/4, 1/8, l/y/6). 

We have numerically studied the RG flows given by 
Eqs. (9) for various starting values of (K,a,b,h). Since 
the Hamiltonian in Eq. (1) does not contain the op- 
erator 3 , we set b = initially. We take a and h 
to be very small initially and sec which set of values 
flows to a non-trivial FP. For instance, starting with 
K = 1, b = 0, and a, h very small, we find that there is 
a line of points which flow to a FP at (K* , a* ,6* , h*) — 
(1.020,0.246,0.122,0.404). This line projected on to the 
(a, h) plane is shown in Fig. 1. We see that K changes 
very little during this flow; if we start with a larger value 
of K initially, then it changes even less as we go to the 
non-trivial FP. It is therefore not a bad approximation 
to ignore the flow of K completely. 
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FIG. 1. RG flow diagram in the (a, h) plane. The solid line 
shows the set of points which flow to the FP at a* — 0.246, 
h* — 0.404 marked by an asterisk. The dotted lines show the 
RG flows in the gapped phases A and B (see text). 

We can characterize the set of points (a, h) lying close 
to the origin which flow to the non-trivial FP. Numeri- 
cally, we find that there is a unique flow line in the (a, h) 
plane for each starting value of K and b — 0, provided 
that a, h are very small initially. This means that a(l) and 
h(l) given by Eqs. (9) must follow the same line regard- 
less of the starting values of a, h. From Eqs. (9), we see 
that if h « a 1 / 2 , then h(l) ~ h(0) exp(2 - K - l/AK)l 
while a(l) <~ cxp(2 — 1/K)l. Hence h must initially scale 
with a as 



h 



(2-K-l/4K)/(2-l/K) 



(11) 



i.e., if K < (1 + V2)/2 = 1.207. For K > 1.207 (i.e., 
A < —0.266), the initial scaling form is given by h ~ a 1 / 2 . 

We now examine the stability of small perturbations 
away from the fixed points. The trivial FP at the 
origin has two unstable directions (a and h), one sta- 
ble direction (b) and one marginal direction (K). The 
non-trivial FP has two stable directions, one unstable 
direction and a marginal direction (which corresponds 
to changing K* and simultaneously a*, b* and h* to 
maintain the relations in Eqs. (10)). The presence 
of two stable directions implies that there is a two- 
dimensional surface of points (in the space of parame- 
ters (a, b, h)) which flows to this FP; the system is gap- 
less on that surface. A perturbation in the unstable di- 
rection produces a gap in the spectrum. For instance, 
at the FP with (K*,a*,b*,h*) = (1,1/4,1/8,1/^6), 
the four RG eigenvalues are given by 1.273 (unstable), 
(marginal), and -1.152 ± 1.067i (both stable). The 
positive eigenvalue corresponds to an unstable direction 
given by (6K,Sa,6b,dh) = <5a(0.113, 1, -0.092, -0.239). 
A small perturbation of size 5a in that direction will 
produce a gap in the spectrum which scales as AE ~ 
|^ a |i/i.273 _ |j a |0.786. correlation length is then given 

by ?~«/AB~ |<H~ ' 786 - 

Figure 1 shows that the set of points which do not flow 
to the non-trivial FP belong to either region A or region 
B. These regions can be reached from the non-trivial FP 
by moving in the unstable direction, with 5a > for 
region A, and 5a < for region B. In region A, the points 
flow to a = oo; this corresponds to a gapped phase in 
which the the xx coupling is larger than the yy and zz 
couplings. In region B, both a and h flow to — oo; this is a 
gapped phase in which the yy coupling is larger than the 
xx and zz couplings. We will now see that the difference 
between these two phases lies in the way in which the 
Z2 symmetry of the Hamiltonian is realized. An order 
parameter which distinguishes between the two phases is 
the staggered magnetization in the y direction, defined 
in terms of a ground state expectation value as 



TO,, 



[ Hm (-1)" < s%sy > ] 1/2 . 



(12) 



This is zero in phase A; hence the Z2 symmetry is unbro- 
ken. In phase B, m y is non-zero, and the Z 2 symmetry 
is broken. The scaling of m y with the perturbation 5a 
can be found as follows [6]. At a — h — 0, the leading 
term in the long-distance equal-time correlation function 
of S y is given by 



(-1)" 
ln n/2K 



(13) 



as we have numerically verified for K — 1. However, Eq. 
(11) is only true if (2 - K - l/4RT)/(2 - 1/K) > 1/2, 



Hence the scaling dimension of S% is 1/4K. In a gapped 
phase in which the correlation length is much larger than 
the lattice spacing, m y will therefore scale with the gap as 
m y ~ (Ai?) 1 / 4 -^ . If we assume that the scaling dimension 
of at the non-trivial FP remains close to 1/4K", then 
the numerical result quoted in the previous paragraph 
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for K = 1 implies that m y ~ |<5a| 0196 for small and 
negative. 

The nature of the transition on the gapless line will 
be discussed in Sec. IV. We will argue there that this 
is a spin-flop transition line. (Spin-flop transitions in 
one-dimensional spin- 1/2 chains have been studied earlier 
[12-14].) 

III. QUANTUM ANNNI MODEL 

We will now apply our results to the one-dimensional 
spin-1/2 quantum ANNNI model [4,7], with nearest 
neighbor ferromagnetic and next-nearest neighbor anti- 
ferromagnetic Ising interactions and a transverse mag- 
netic field. The Hamiltonian is given by 

H A = £ [ - 2J 1 T„*T n x +i + + ] , 

n 

(14) 

where Ji, J2 > 0, and the T" are spin-1/2 operators; 
we can assume without loss of generality that T > 0. 
The quantum Hamiltonian in Eq. (14) is related to the 
transfer matrix of the two-dimensional classical ANNNI 
model; the finite temperature critical points of the latter 
are related to the ground state quantum critical points 
of Eq. (14), with the temperature T being related to the 
magnetic field T. 

Some earlier studies showed that the model has a float- 
ing phase of finite width which is gapless [4]. A re- 
cent bosonization study reached the same conclusions 
[7]. (Recent numerical studies of the two-dimensional 
classical ANNNI model at finite temperature have led to 
contradictory results for the width of the floating phase 
[15].) All these studies (both analytical and numerical) 
indicate that the phase transition is of the Kostcrlitz- 
Thouless type (with £ diverging exponentially) from the 
high-temperature side (i.e., from region B in Fig. 1 for 
the quantum ANNNI model), and is of the Pokrovsky- 
Talapov type [16] (with £ diverging as a power-law) from 
the low-temperature side (i.e., from region A in Fig. 1). 

We will now apply our results to the quantum ANNNI 
model. Consider a Hamiltonian which is dual to Eq. (14) 
for spin-1/2; this will turn out to be a special case of our 
earlier model. The dual Hamiltonian is given by [4,17] 

H D =Y,\ J2S x n S x n+1 + rS« n S v n+1 - JxS* ] , (15) 

n 

where S 1 " are the spin-1/2 operators dual to T" (for in- 
stance, SI = 2T*T* +1 andT r f = 2S*£_ 1 S£). After scaling 
this Hamiltonian by an appropriate factor, we see that it 
has the same form as in Eq. (1), with 

J2-T 

a = jT+t ' 



A = . (16) 

Hence it follows that the quantum ANNNI model has a 
line of points in the (J2/J1, T / J\) plane on which the 
system is gapless. From Eq. (11), we see that the shape 
of this line is given by J\ ~ ( J 2 — T) 3 / 4 as J\ — > 0. 

The analysis in Sec. II indicates that as the transition 
line is approached, £ should diverge as a power-law from 
both sides. We now have to reconcile this with some 
of the earlier analytical [4,7] and numerical [15] studies 
which showed that as one approaches the floating phase, 
£ diverges as a power-law from phase B but exponentially 
from phase A. The important point is that these earlier 
studies were carried out at values of J2 / Ji which are close 
to 1, while our RG results are expected to be valid only 
if a, h are small, i.e., if J2/J1 is large. If J2/J1 is close to 
1 , the situation is quite different for the following reason. 
Exactly at J2/J1 = 1 and T = 0, the Hamiltonian in Eq. 
(15) can be written in the form 

H M c = -h ^2 ( S n ~ 2-* ~~ 2-* ' ^ 

n 

This is a multicritical point with a ground state degener- 
acy growing exponentially with the system size, since any 
state in which every pair of neighboring sites (n, n+1) has 
at least one site with S x = 1/2 is a ground state. We can 
now study what happens when we go slightly away from 
this multicritical point. To lowest order, this involves 
doing perturbation theory within the large space of de- 
generate states. An argument due to Villain and Bak [4] 
shows that if J 2 — J\ and Y are non-zero but small, then 
the low-energy properties of Eq. (15) do not change if 
TSlS v n+l is replaced by (r/2)(S*Sj[ +1 + S*S* n+1 ). (This 
is because the difference between the two kinds of terms 
is given by operators which, acting on one of the degen- 
erate ground states, take it out of the degenerate space 
to a higher excited state in which a pair of neighboring 
sites have S x = —1/2.) Thus the fully anisotropic model 
becomes equivalent to a different model which is invari- 
ant under the U(l ) symmetry of rotations in the y — z 
plane. The U(l) symmetric model has been studied ear- 
lier using bosonization [2,11,18]; it has a gapless phase of 
finite width which lies between two gapped phases. Thus 
the difference between our study (in which J 2 — T and J\ 
are small) and the earlier studies (in which J 2 — J\ and 
r are small) is that they have different symmetries away 
from the transition line, namely, Z 2 and U(l) respec- 
tively. Our study and the earlier studies are therefore 
complimentary to each other; a combination of the two 
leads to a complete understanding of the model over the 
entire parameter range. 

To summarize, the transition from phase A to phase B 
can occur either through a gapless line (if o, h are small) , 
or through a gapless phase of finite width (if a, h are 
large). The complete phase diagram of the ground state 
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of Eq. (15) is shown in Fig. 2 [4]. The three major 
phases shown are distinguished by the following prop- 
erties of the expectation values of the different compo- 
nents of the spins. In the antifcrromagnctic phase, the 
spins point alternately along the x and —x directions. 
In the spin- flop phase, they point alternately along the 
y and — y directions, with an uniform tilt towards the x 
direction. In the ferromagnetic phase, all the spins point 
predominantly in the x direction. The antiferromagnetic 
and spin-flop phases are separated by a floating phase 
of finite width for J2/J1 close to 1/2, and by a spin-flop 
transition line for large values of J 2 j J\. We conjecture 
that the floating phase and the spin-flop transition line 
are separated by a Lifshitz point as indicated in Fig. 2. 
The disorder line and the Ising transition are discussed 
in Sees. V and VI respectively. 

We should point out here that in terms of the original 
Hamiltonian in Eq. (14), some of the phases shown in 
Fig. 2 have somewhat different names [4]. The spin- flop 
phase is called the paramagnetic phase; this is further 
divided into two phases by the disorder line, namely, a 
commensurate phase to the left and an incommensurate 
phase on the right of the disorder line. The antiferro- 
magnetic phase is called the antiphase. 
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FIG. 2. Schematic phase diagram of the model described 
in Eq. (15). The various phases and transition lines are ex- 
plained in the text. The initials FP, KT and PT stand for 
floating phase, Kosterlitz-Thouless and Pokrovsky-Talapov 
respectively. 



IV. CLASSICAL LIMIT 

In this section, we would like to provide a physical 
picture of the gapless line in the (a, h) plane by looking at 
the classical limit of Eq. (1). Consider the Hamiltonian 

Hsi = £ [(l + a)S*S* n+1 + (l-a)« +1 

n 

+ A - 2Sh SI } , (18) 

where the spins satisfy = S(S + f ), and we are in- 
terested in the classical limit S — > 00 [17]. (We have 
multiplied the magnetic field by a factor of 2S in Eq. 
(18) so that we recover Eq. (1) for spin-1/2.) We assume 
as before that the zz coupling is smaller in magnitude 
than the yy coupling. Then the classical ground state 
of Eq. (18) is given by a configuration in which all the 
spins lie in the x — y plane, with the spins on odd and 
even numbered sites pointing respectively at an angle of 
a\ and — a 2 with respect to the x-axis. The ground state 
energy per site is 

e(ai,a 2 ) = S 2 [— h (cos ai + cos a 2 ) + cos(ai + 0:2) 
+ a cos(ai — a 2 )]. (19) 

Minimizing this with respect to ot\ and a 2 , we discover 
that there is a special line given by h 2 = 4a on which 
all solutions of the equation 

ai-a 2 . a 1 +a 2 , . 

h cos( ) = 2 cos( ) , (20) 

give the same ground state energy per site, e = — (1 + 
a)S 2 . The solutions of Eq. (20) range from ot\ = a 2 = 
cos _1 (/i/2) to ct\ — ir,a 2 = (or vice versa); in the 
ground state phase diagram of the ANNNI model, these 
two configurations correspond respectively to a antiferro- 
magnetic alignment of the spins with respect to the y-axis 
(with a small tilt towards the x-axis if h is small), and 
an antiferromagnetic alignment of the spins with respect 
to the x-axis. The curve h 2 = 4a is therefore a phase 
transition line, and the form of the ground states on the 
two sides shows that there is a spin-flop transition across 
that line. Further, we see that for h 2 = 4a, there is a one- 
parameter set of classical ground states (characterized by, 
say, the value of ct\ which can go all the way from to 
2-7T in the solutions of Eq. (20)) which are all degenerate. 
Hence the symmetry is enhanced from a Z 2 symmetry 
away from the line to a U(l) symmetry (of rotations in 
the x — y plane) on the line. We therefore expect a gap- 
less mode in the excitation spectrum corresponding to 
the Goldstone mode of the broken continuous symmetry. 
We can find this gapless mode explicitly by going to the 
next order in a 1/S expansion [17]. 

The above arguments provide some understanding of 
why one may expect such a gapless line in the spin-1/2 
model also. Note however that the bosonization anal- 
ysis gives the scaling form in Eq. (ff) for h versus a; 
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this agrees with the classical form only if A < —0.266. 
Further, in the classical limit, the transition across the 
gapless line is of first order, whereas it is of second order 
in the spin- 1/2 case. There is probably a critical value 
of the spin 5 above which the transition is of first order. 
(For the U(l) symmetric model described by Eq. (21) 
below, it is known that the transition is of first order if 
5> 1 [14].) 

The classical limit also makes it clear why our model 
has a different behavior from the U(l) symmetric model 
governed by the Hamiltonian 

H S 2 = £ [(l + a)S:s: +1 + (l-a)SlSl +1 

n 

+ (l-a)S*5* +1 - 2ShS x n \. (21) 

In the limit 5 — > oo, there is now a two-parameter set 
of degenerate ground states on the line h 2 = 4a; these 
are obtained by taking the one-parameter family of con- 
figurations given in Eq. (19) and rotating them by an 
arbitrary angle about the x-axis. Hence, the symmetry 
of this model is enhanced from U (1) to SU (2) on the line 
h 2 = 4a, and there are now two Goldstone modes in- 
stead of one. Considering this difference in symmetry for 
large 5, it is not surprising that even the spin- 1/2 models 
with 1/(1) symmetry and Z 2 symmetry respectively ex- 
hibit very different behaviors at the spin-flop transition 
line. 

V. DISORDER LINE 

We have seen that as the magnetic field h is increased 
from zero for the spin-1/2 model described by Eq. (1), 
there is a spin-flop transition at a critical field h c whose 
value depends on a and A. One might wonder what 
happens if the field is increased well beyond h c . 

It turns out that above h c , there is an interesting value 
of the field h = hd where the ground state of the model 
is exactly solvable [19,20]. This field is given by 

h d = v/2(l + a + A) . (22) 

At this point, the ground state has a very simple direct 
product form in which all the spins lie in the x — y plane, 
with the spins on even and odd sublattices pointing at the 
angles a and —a respectively with respect to the x-axis, 
where 

a = cos" 1 ^) . (23) 

To show that this configuration is the ground state of 
the Hamiltonian, we observe that the Hamiltonian can 
be written, up to a constant, as the following sum 

,2n+l + #2n,2n-l] ; 

n 



H 2n ,2n±i = (cos a S 2n + sin a 5f„ - -) 

x(cosa 5f n±1 + sin a S y 2n±1 - ^) 
+ (cosa 5f„-sina S* n - 1) 
x(cosa 5f„ ±1 - sin a S v 2n±1 - 1) 

+ A [ 1 - (cos a 5 2 *„ + sin a S v 2n ) 

x (cos a 5f„ ±1 - sin a S% n±1 ) 
- (sin a 5f„ - cos a 5fJ 
x (sin a 5f„ ±1 + cosa S v 2n±1 ) 
+ SLSI n±1 } , (24) 

where a is given in Eqs. (22-23) . We now use the theorem 
that the ground state energy of H is greater than or equal 
to the sum of the ground state energies of H 2n ,2n±i, with 
equality holding if and only if there is a state which is 
simultaneously an eigenstate of all the H 2n ^ 2n ±\. Now, 
each of the Hamiltonians H 2n _ 2n ±\ in Eq. (24) is a sum 
of three operators whose eigenvalues are non-negative if 
A > [20]. The state described in Eq. (23), in which 
all the spins on the even sublattice satisfy cosa5f„ + 
sin aSf^ = 1/2 and all the spins on the odd sublattice 
satisfy cos aS 2n+1 — sin aS 2n+1 = 1/2, is the ground state 
of all the Hamiltonians in Eq. (24) with zero eigenvalue. 
We can actually show, by looking at a two-site system 
governed by a single Hamiltonian H 2rii2n+ i, that even 
if A < 0, the state described above is its ground state 
provided that 1 — a > —A, i.e., as long as the magnitude 
of the zz coupling is smaller than the yy, which is what 
we have assumed already. 

For a given value of A, the line in the (a, h) plane 
described by Eq. (22) is called a disorder line because 
the direct product form of the ground state implies that 
the two-spin correlation function < S^S^ > — < S" >< 
Sm > (with a,(3 = x, y, z) is exactly zero \im^n. Hence 
the correlation length is extremely short. The disorder 
line exists even for values of the spin larger than 1/2. 
Starting with the Hamiltonian in Eq. (18), one finds a 
disorder line at the same value of h given in Eq. (22). 
The proof that it is a disorder line is similar to the proof 
given above for the spin-1/2 case if A > 0. We will not 
study here how far the proof can be extended to negative 
values of A; for spin S, this requires an examination of 
the spectrum of a two-site problem governed by a (25 + 
1) x (25 + 1) dimensional Hamiltonian matrix. 

VI. ISING TRANSITION 

If the magnetic field h is increased even further, the 
system undergoes an Ising transition [4]. If the yy and 
zz couplings are equal (i.e., 1 — a = A), this occurs at 
a saturation field h s = 2, where there is transition to a 



6 



state in which all the spins point along the x-axis. But if 
the yy and zz couplings are not equal, there is no satura- 
tion of the spins for any finite value of the field although 
the ground state expectation value of S% approaches 1/2 
(as (1 — a — A) 2 /h 2 ) as h goes to infinity. (This can be 
shown by considering a two-site system and doing per- 
turbation theory in the limit h — ► oo.) However, there 
is still a transition field h s beyond which a Z 2 symmetry 
of a different kind is broken. To see this, we consider a 
Hamiltonian H which is dual to the Hamiltonian given 
in Eq. (1). This is given by 

H = J2 [{l + a)T*T* +2 + 

n 

- 2A T*_{FITI +1 - 2h T:t: +1 } . (25) 

This Hamiltonian is invariant under the global Z 2 trans- 
formation T% - —T%, T% - T%, T* - -T*. For A - 0, 
this Z 2 symmetry is known to be broken if h is larger 
than a critical value h s [4]. We expect that this is will 
be true even if A ^ 0. The order parameter for this 
symmetry is 

m x = [ lim < T*T: > ] 1/2 . (26) 

n — >oo 

Note that in terms of the operators TqT% is equal 
to a string of operators, (1/4) I]m=o(^m)' Similarly, the 
order parameter (— 1)"5q5^ in Eq. (12) is equal to the 
string of operators ((-l)»/4) FEUiM)- 

VII. DISCUSSION 

We have shown in this paper that the XYZ spin- 1/2 
chain in a magnetic field exhibits a gapless phase on a 
particular line. It would be interesting to use numerical 
techniques like the density-matrix renormalization group 
method [21] to examine various ground state properties 
of this model, in particular, to study the behavior of the 
order parameter defined in Eq. (12), and to find out if 
there is indeed a Lifshitz point as conjectured in Fig. 2. 

Finally, the RG equations studied in this paper appear 
in other strongly correlated systems, such as the prob- 
lem of two spinlcss Tomonaga-Luttinger chains with both 
one- and two-particle interchain hoppings [10], and one- 
dimensional conductors with spin-anisotropic electron in- 
teractions [11]. The gapless phase may therefore also ap- 
pear in other systems. 
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